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Abstract 

Lie-algebraic and quantum-algebraic techniques are used in the analysis of ther- 
modynamic properties of molecules and solids. The local anharmonic effects are 
described by a Morse-like potential associated with the su(2) algebra. A vibrational 
high-temperature partition function and the related thermodynamic potentials are 
derived in terms of the parameters of the model. Quantum analogues of anhar- 
monic bosons, g-bosons, are introduced and used to describe anharmonic properties 
of molecules and solids. It is shown that the quantum deformation parameter is 
related to the fixed number of anharmonic bosons and the shape of the anharmonic 
potential. A new algebraic realization of the (/-bosons, for the case of q being a root 
of unity is given. This realization represents the symmetry of a linear lattice with 
periodic boundary conditions. 

1 Introduction 

The algebraic approach has been developed as an alternative to ah initio and Dunham-like 
approaches to describe the molecular vibrational degrees of freedom Jl|, [2| . The algebraic 
approach provides the energy as an analytic function of the quantum numbers. The 
Hamiltonian is written as an algebraic operator using the appropriate Lie algebras. The 
technical advantage of the algebraic approach is the comparative ease of the algebraic oper- 
ations. Equally important is the result obtained by comparison with the experiment, that 
there are generic forms of algebraic Hamiltonians and that entire class of molecules can be 
described by common Hamiltonian where only the parameters are different for different 
molecules. In its initial stage of development [@]-@, the algebraic approach has sought 
to show why and how it provides a framework for the understanding of large-amplitude 
anharmonic motion. The anharmonicities are introduced by means of dynamical groups 
that correspond to anharmonic potentials and which constrain the total number of levels 
to a finite value. The current algebraic models ff|-f| combine Lie algebraic techniques, de- 
scribing the interatomic interactions, with discrete symmetry techniques associated with 
the local symmetry of the molecules. In the anharmonic oscillator symmetry model 0, 
the local internal coordinates are given in terms of u(2) algebras. The u(2) interactions 
correspond to anharmonic coupling of anharmonic oscillators which approximate the in- 
teractions between the Morse oscillators. 
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The algebraic anharmonic model has been developed to analyze molecular vibrational 
spectra 0-0. It provides a systematic procedure for studying vibrational excitations in 
a simple form by describing the stretching and bending modes in a unified scheme based 
on u(2) algebras which incorporate the anharmonicity at the local level. 

The success of the algebraic models in the analysis of molecular vibrational spectra 
has led to the development of similar models for molecular thermodynamics. Kusnetzov 



10 has used an algebraic approach, based on approximation of the classical density of 



states, to study thermodynamic properties of polyatomic molecules. Angelova and Frank 
|TTi |T2| , p~3|] have applied the algebraic model |], to the vibrational high-temperature 



thermodynamics of diatomic molecules and derived the vibrational partition function and 
the important thermodynamic functions, such as mean energy and specific heat, in terms 
of the parameters of the model. 

The Morse-like potential, which represents the anharmonicities at the local level, leads 



to a deformation of the harmonic oscillator algebra. Angelova, Dobrev and Frank JT5| 
have derived this deformation using quantum analogue of the anharmonic oscillator. We 
have described the anharmonic vibrations as anharmonic g-bosons using first-order of the 
expansion of a quantum deformation and found relations between the parameters of the 
algebraic model and the quantum deformation parameter. 

The aim of this paper is to review the applications of the classical and quantum 
algebras to molecular thermodynamics. In Section 2, the framework of the algebraic model 
is given. The vibrational partition function and the related thermodynamic functions, 
such as mean energy, specific heat and the mean number of the anharmonic bosons, are 
discussed in Section 3. The anharmonic g-bosons are discussed in Section 4. A g-bosonic 
deformation of first order is considered and it is shown that the corresponding quantum 
deformation parameter is related to the shape of the anharmonic potential well and the 
fixed number of anharmonic bosons. The g-bosons at roots of unity, which give rise to 
a finite-dimensional periodic structure are discussed and their applications to a linear 
lattice with periodic boundary conditions are given. 



2 Algebraic Model 

The algebraic model exploits the isomorphism of su(2) algebra and the one-dimensional 
Morse oscillator. The one-dimensional Morse Hamiltonian can be written in terms of the 
generators of su(2), 

H M = j (A" 2 - 4Jf) = ^{Jj- +J-J+- AO (1) 

where A is a constant dependent on the parameters of the Morse potential. The eigen- 

states, |[A/],u), 



mv) = \i^ 7 ^(j-) v \w\,o) (2) 
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correspond to the u(2) D su(2) symmetry-adapted basis, where M is the total number 
of bosons fixed by the potential shape, and v is the number of quanta in the oscillator, 
v = 1, 2, . . . , [y] . These wave functions can be further symmetry-adapted to the local 
symetry of the molecules in a way described for example in |7| . 

The value of M is dependent on the depth D and the width d of the Morse potential 
well 0, 



8fiDd 2 



(3) 



where /i is the mass of the oscillator. The parameters A and M are related to the usual 
harmonic and anharmonic constants uj e and x e u e used in spectroscopy f|, H| , 



A * 



2D 



(4) 
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The anharmonic effects are described by anharmonic boson operators [lj] 

? J+ t* J- 



Jf 



2 



J, 



(5) 



where v is the Morse phonon operator with an eigenvalue v. The operators satisfy the 
commutation relations, 



6, v 



b, b 



2v 



oo, in which case 



b,V 



The harmonic limit is obtained when TV 
usual boson commutation relations. 

The Morse Hamiltonian can be written in terms of the operators b and b\ 



(6) 

1 giving the 



Hm ~ 

which corresponds to vibrational energies 



x - (IV + Vb 



1 v 2 



1,2, 



J£ 
2 



(7) 



(8) 



where uq is the harmonic oscillator frequency. Thus, the spectrum of the Morse potential 
leads to a deformation of the harmonic oscillator algebra. A more detailed relationship 
between the Morse coordinates and momenta and the su(2) generators can be derived 
through a comparison of their matrix elements 0. Note that for an infinite potential 
depth, Af— >oo, the Morse potential cannot be distinguished from the harmonic potential. 
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3 Thermodynamic Vibrational Functions 



3.1 Vibrational Partition Function 

In the anharmonic algebraic approach, the vibrational partition function of a diatomic 
anharmonic molecule is 

[A/72] 

Z M = E (9) 

v=0 

where (3 = 1/ksT, the vibrational energies e v are given by equation (|j) and M is the fixed 
total number of anharmonic bosons discussed in the previous section. Introducing new 
parameters, a = ^p-, A/o = [y] and / = Af Q — v, the exact value of vibrational partition 
function can be written as, 

A/o 

Z A r = e- a ^e"^( A '' , ^ /2 ). (10) 
1=0 

At high temperatures T, for Ao large and a small, the sum can be replaced by the integral, 

VaMo 

W— e-^ Vo+ij / e s "ds (II; 




where s = w -^Z. This integral can be evaluated exactly in terms of the error function, 
erf i (-y/oA/o) (as defined in fll6|), 

^ = iy / ^e^^ +1 )erfz( v ^) . (12) 

Equation ([T2] ) represents the high-temperature value of the vibrational partition function 
in the Morse-like spectrum. The partition function is expressed in terms of the parameters 
of the algebraic model Ao and a. When A/q — > 00, the harmonic limit of the model is 
obtained, 

Af e~ a e~ a T _^ 

Zoo ~ 77 ~ = T^e 2T (13) 

2aAT -l 2a 9 v ; 

which coincides with the harmonic vibrational partition function of a diatomic molecule at 
high temperatures. Here, = hooo/kB, is the usual characteristic vibrational temperature 
of the molecule and the model parameter a = G/2T. The expression for the partition 
function ( |T2"D can be further generalised for polyatomic molecules. 



Using the values of the harmonic and anharmonic constants MM for the zero lines of 



the diatomic molecule 1 H 35 C1 and the equations (H) and @), we obtain the total number 
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Figure 1: Vibrational partition function Z 56 as a function of a. 



of anharmonic bosons, fixed by the shape of the Morse potential, M = 56, and the total 
number of quanta in the oscillator, A/"o = 28. The characteristic vibrational temperature 
of the molecule is = 4300K. 



Substituting the value of A/q = 28 in equation ( |12"1) we can calculate the partition 
function, Z 56 , for the molecule 1 H 35 C1 as a function of the parameter a. The anharmonic 
effects are essential at the high temperatures T > Q, i.e. a < 0.5. The graph on Figure [l| 
represents the partition function Z 56 given by equation fll^D for the values of the parameter 
a between and 0.5 (solid line). The exact partition function from equation (pTDj) is given 
for comparison (dashed line). It is clear, that the integral approximation is in a very good 
agreement with the exact representation and does not change the value and appearance 
of the partition function. 



3.2 Mean Vibrational Energy 

The mean vibrational energy is given by 



U N = -^\nZ N 



Taking into account that 



2Z N da 



dZ M 



ZM -(N- + l)Z J , + Moe ~ 



(14) 



(15) 



da 2a " ' 2a 

we obtain the following expression for the mean vibrational energy in terms of the partition 
function Zyv, 



1 AToe- a \ 
2a 2aZj^f ) 



(16) 
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Substituting Z^ with equation (0) gives the following expression for the mean energy, 
U/s, in terms of the parameter a, 

UM = ^( l+ „ t+ iM_jfLJ\. (17) 

2 I 2a V a7rerf« (VaA/o) y 

The harmonic limit is obtained from equation (0), when Ao — > oo and Zjv is given by 

f/oo ~ — (1 + -) ~ — + fcsT. (18) 
This is the classical mean energy of a diatomic molecule at very high temperatures. 

3.3 Specific Heat 

The vibrational part of the specific heat is, 

_ dU N huj dU N 
CN -^T-~2k^^a~- (19) 

Substituting Uj^ with equation (|T6|) and using equation QT5T) , we obtain 

The equation (|20| ) represents the vibrational specific heat in the algebraic model in terms 
of the partition function Z^. Substituting Z^ with (|P2"|) in the equation (p0|) , we obtain 
the dependence of the specific heat CV on the parameters of the model a and Ao, 



T *"V ~ er fi P° - 2 - V — erf i J (21) 

When Ao — > oo, the harmonic limit of the model gives the vibrational specific heat 
of a diatomic molecule at very high temperatures, 

Oc ~ k B . (22) 
3.4 Mean Number of Anharmonic Bosons 

The mean vibrational energy in the anharmonic model can be written in terms of mean 
number of anharmonic quanta, each with energy huo, 

U N = hu ((ujsf) + i J (23) 
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Figure 2: Mean Number of Anharmonic Bosons z/ 56 as a function of a 
Substituting Uj^ by equation fll6|) , we obtain {ujf) in terms of the partition function Zjy, 



M = — + 



oe 



2 4a AaZtf 

Using the expression (|l|) in equation (f24"D, we obtain the high-temperature value, 



(24) 



= IT + T~ 



3 aA/b 



2 4a V 47ra erf z ( v / aA/" .) 



(25) 



The harmonic limit is obtained from equation fl24|) when A/o — > oo and zT^ is given by 
the expression (|13|), 



(foo) 



(26) 



The graph of the function (z/56) for the molecule 1 H 35 C1 is given on Figure |2| (solid line). 
The harmonic limit (voo) is given for comparison (dashed line). 



4 Applications of q-Bosons 
4.1 Anharmonic g-Bosons 

We have shown in \TE\ that the anharmonic bosons b, from equations (|5|) can be obtained 



as an approximation of the g-bosons [17], [18|, [19|. The g-bosons, enter the Heisenberg-Weyl 
g-algebra HW q by the following commutation relations: 



(27) 
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where the deformation parameter q is in general a complex number. When q — 1, the 
boson commutation relations of the harmonic oscillator are recovered. 
A possible Hamiltonian for the system (^7]) is: 

H = -(aat + a t a ) = -C + ^ -. (28) 

2 V ; 2 2 q-1 K ' 

where the Casimir operator C can be written in the form: 



t t + g" - 2 
C = aa> + a'a 



(29) 

q- 1 

and satisfies the commutation relations, 

[C, a] = [C, a f ] = [C, n] = (30) 



As shown in [15|, the anharmonic bosons ((g) can be obtained from the g-bosons 
for real values of the deformation q close to 1, q < 1, using an expansion of first order of 
q in terms of a parameter p, p = 1/(1 — q), 

q h = l-- . (31) 
p 

where l/p<Cl. If we now substitute the approximation for q n from equation (|3~I| ) in the 



commutation relations (|27 ) and identify the parameter p with Af/2, n with ?) and the 
creation and annihilation operators a,a\ with b,tf respectively, we recover the su{2) 
anharmonic relations (§). 

For q < 1, the case of harmonic and anharmonic vibrations in molecules and solids 
is retrieved. The form (|]) of the sw(2) commutation relations can be considered as a 
deformation of the usual (harmonic oscillator) commutation relations, with a quantum 
deformation parameter p = J\f/2. 

This gives a physical interpretation of the quantum deformation. The quantum defor- 
mation parameter p is the fixed number Afo of the anharmonic bosons in the oscillator. 
Using the relation between the fixed number of anharmonic bosons Af and the characteris- 
tics of the Morse potential ([5|), we can conclude that the quantum deformation parameter 
is also determined by the depth, the width and in general the shape of the Morse poten- 
tial well. For the molecule ^Cl, p = 28 and q = 27/28. Substituting jV = p in the 
expressions (0), Q17D , (pT|)and fl25|), we obtain the thermodynamic properties of diatomic 
molecules in terms of the deformation parameter p. 

The case q > 1 is also very interesting and is related to Bose-Einstein condensation 



and superfluidity 24, 25 



4.2 q-Bosons at roots of unity 

In ||15|| , we have presented a new algebraic realization of the g-bosons , for the case of q 
being a root of unity, q N = 1 (integer N > 1), which corresponds to a periodic structure 
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described by a finite-dimensional representation. Note that the integer N in this section 
is not related in general to the total number of bosons M in the previous sections. 
We consider the operator K = q n with commutation relations, 



[a, a 



r l-K, Ka = q~ 1 aK, Ka) = qa ] K. (32) 



| ) is the vacuum state which is annihilated by the operators lowering the boson number 
and is an eigenvector of the number operator: 

a|0) =0 , K\) = <f|0), (or n|) = /x|0» (33) 

where \x in the generic case is an arbitrary complex number. The states of the system are 
built by applying the operators raising the boson number: 

\k) = (at) fc |0). (34) 

The action of the algebra on the basis \k) is: 



K\k) 


= q" +k \k) 


a\k) : 


Q k — 1 


q-1 


a*\k) 


= \k + l). 



(35) 

The representation space, V^, is infinite-dimensional for a generic deformation parameter. 

Now, let q is a nontrivial root of unity, i.e., q — 1 for integer N > 1. In this case we 
have: 

q N - 1 

a \N) = q» ^— — \N - 1) = 0. (36) 

All states with k > N form an infinite dimensional invariant subspace, 1^. We then 
obtain a finite-dimensional representation space, which is the factor-space, = V^/I^, 
with dimension N. Considering the action of the algebra on the states, \k + mN), for 
fixed k < N and for all non-negative integer m, we have shown in that the structure 
is periodic, i.e., the action of the algebra on all states \k + mN) (for fixed k) coincides. 
Thus, it is sufficient to consider the states \k) with k < N . Let us denote these identified 
states by \k), k = 0, . . . , N — 1. They form a finite-dimensional representation space F M 
with dimension N. The action of the algebra on these states is: 

K\k) =q fl \k) 

ojjfe) =q^^—^\k^l) 

Ss 1 

a f = \k + 1) , fc<iV-l 

^ |iV~l) = |0), |0) = |iV). (37) 
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This represents a finite-dimensional system, where the boson number lowering operator 
acts in the usual way, (in particular, it annihilates the vacuum state |0)), but the boson 
raising operator acts cyclicly. It has a non-zero action on all states and the vacuum state 
may be obtained not only by the action of the lowering operator but also by the action 
of the raising operator producing a jump from \N — 1) to |0). One realization of this 
operator is a two-level system, obtained for N = 2 (equivalent to q = — 1, fermions). For 
N > 2, systems with finite number of levels and population inversion are illustrations of 
possible action of these operators. In the objects "quons" (q N = 1) are introduced as 
interpolating between fermions (N = 2) and bosons(iV — > oo). 

For large N periodicity of the type decribed above appears in crystals. We have 
shown that it represents the periodic boundary conditions, first proposed by Born and 
von Karman ||27]| . The periodic boundary conditions |28| are imposed on the translational 



symmetry, which strictly speaking is a property of an infinite crystalline lattice, to allow 
its use for finite crystals. The periodic boundary conditions determine the number of 
the allowed wave- vector states in the Brillouin zone model and imply additional selection 
rules on certain frequencies. 

For the classic example of a linear lattice of identical particles with periodic bound- 
ary conditions, the equilibrium positions of the particles are given by t n = nt, n = 
0, 1, . . . ,N — 1, and the periodic boundary condition requires tjv = t = 0., where t is the 
vector of primitive translations and N is a large positive number. 

The symmetry operations of the linear lattice form a cyclic finite group of order N 
with a generator, the primitive translation {-E^t}. Here, the Seitz notation is used to 
represent a translation and E is the identity, E = {E\0}. 

The n-th element of the group is 

{E\t} n = {E\t n }, n = l,2,... ,N-1. (38) 
The product of two elements of the group is an element of the group 

{E\t m }{E\t n } = {E\t m+n }, m, n = 1, 2, . . . , N - 1, (39) 
where m + n = (m + n) (mod N). The identity is 

{E\t} N = {E\t } = {E\0}. (40) 

Using the action ( j37|) of the operator a* on the states |fc), we have shown that the raising 
operator a) is isomorphic to the generator {-E^t}. 

Thus, the symmetry group of the lattice with periodic boundary conditions is isomor- 
phic to the finite cyclic group of order iV with a generator, the operator a'. This group 
can be used with the other symmetry operations of one-dimensional crystalline or polymer 
Hamiltonians. The boundary conditions can be generalised for the three-dimensional case 
by introducing raising operators for each dimension. 
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5 Conclusion 



In this paper, we have used the classical and quantum algebras to study the vibrational 
thermodynamic properties of molecules and solids. The algebraic approach is applied 
to those thermodynamic properties of diatomic molecules which at high temperatures 
strongly depend on the anharmonic effects. The vibrational thermodynamic functions, 
such as partition function, mean energy and specific heat, are derived in terms of the 
parameters of the algebraic model and their properties are discussed. 

The application of a g-algebra to physical problems is often lacking an appropriate 
interpretation for the deformation parameters and often applications are carried out where 
generalisation to g-deformed versions of well known models are made with no simple 
interpretation. We have shown that the 1/p approximation leads to the su(2) algebra 
and to an interpretation of p in terms of the Morse potential anharmonicity. We have 
found a physical interpretation of the quantum deformation, showing that the deformation 
parameter is related to the fixed number of anharmonic bosons and the shape of the Morse 
potential. We have found that all vibrational thermodynamic properties (e.g. mean 
energy, specific heat, mean number of anharmonic bosons) of the molecules depend on 
the corresponding quantum deformation parameter. 

A new application of the HW q algebra when q is a root of unity is discussed, which 
gives a periodic structure described by a finite-dimensional representation. The raising 
operator belonging to this structure generates a group isomorphic to the symmetry group 
of a linear lattice with periodic boundary conditions. The latter may provide a useful 
framework for the deformation of crystalline or polymer Hamiltonians. 
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